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ABSTRACT 

The Sunyaev-Zel'dovich (SZ) effect is a direct probe of thermal energy content of the 
Universe, induced in the cosmic microwave background (CMB) sky through scattering 
of CMB photons off hot electrons in the intracluster medium (ICM). We report a 9<r 
detection of the SZ signal in the CMB maps of Wilkinson Microwave Anisotropy Probe 
(WMAP) 3yr data, through study of a sample of 193 massive galaxy clusters with 
observed X-ray temperatures greater than 3 keV. For the first time, we make a model- 
independent measurement of the pressure profile in the outskirts of the ICM, and 
show that it closely follows the profiles obtained by X-ray observations and numerical 
simulations. Wc find that our measurements of the SZ effect would account for only 
half of the thermal energy of the cluster, if all the cluster baryons were in the hot ICM 
phase. Our measurements indicate that a significant fraction, 35% ± 8%, of baryonic 
mass are missing from the hot ICM, and thus must have cooled to form galaxies, 
intracluster stars, or an unknown cold phase of the ICM. There does not seem to be 
enough mass in the form of stars or cold gas in the cluster galaxies or intracluster 
space, signaling the need for a yet-unknown baryonic component (at 3er level), or 
otherwise new astrophysical processes in the ICM. 



1 INTRODUCTION 



Galaxy clusters are remarkable laboratories for studying 
structure formation and cosmology. Thanks to their physi- 
cal size of ~ Mpc, they can be easily resolved in different 
frequencies, ranging from radio, optical to X-rays, and thus 
any model for cluster physics can be independently tested 
against various observations. The simple model that is borne 
out of these studies can then be used as a standard can- 
dle/ruler that probes the geometry an d dynamics of the 
Universe at cosmological redshifts (e.g., lBorganill2006l ). 

However, the simple theoretical picture of cluster forma- 
tion through gravitational collapse is complicated by other 
astrophysical processes such as shock-heating and radiative 
cooling of gas, star formation, energy feedback and chem- 
ical enrichment of the ICM by supernovae and accreting 
super-massive black holes. Most of our observational under- 
standing of the intracluster medium (ICM) to date has come 
from the study of its diffuse X-ray emission, which has led 
to a multitude of theoretical investigations on the eff ects 
of gas cooling and heating of the ICM (see IVoitl 120051 . for 
a review). Although modern high-resolution hydrodynamic 
simulations have enabled detailed theoretical modeling of 
cluster formation, the details of star formation or interac- 
tions with super-massive black holes are still far beyond the 
numerical resolution limit, and our understanding of details 
and relative importance of most of these processes still re- 
main quite uncertain. Therefore, further observational stud- 



ies of clusters in X-ray and other wavebands will continue to 
be the key in guiding physical modeling of cluster properties 
and formation. 

The thermal Sunyaev Zel'dovich (SZ) effect is a unique 
observational probe of galaxy clusters detectable in ra- 
dio and microwave frequencies. It is a spectral distor- 
tion in CMB spectrum caused by inverse-Compton scat- 
tering of CMB photons off hot electrons in the ICM 
l|Sunvaev fc Zel'dovichl 19721 ; iBirkinsha^ll 1999ft . As the num- 
ber of photons does not change through scattering, the SZ ef- 
fect causes a decrement (increment) in the intensity of CMB 
at low (high) frequencies. The effect has the unique prop- 
erty that its signal is independent of redshift, and it directly 
probes the ICM pressure, or thermal energy density. These 
features make the SZ effect a particularly unique and pow- 
erful observational tool for detecting clusters at cosmologi- 
cal redshifts [z > 0.5) and hence a promising cos mological 
probe of dark energy (see ICarlstrom et alj |2002| . for a re- 
view) . 

Although various CMB experiments are now aiming 
at carrying out SZ cluster surveys, all-sky CMB surveys 
with sufficient resoluti on to detect many SZ clu sters di- 
rectly (such as Planck; [Ccisbucsch & Hobson 20(3) are yet 
non-existent. The highest resolution all-sky map of the 
CMB sky has been recently released by the Wilkinson Mi- 
crow ave Anisotropy Probe (WMAP) group (jHinshaw et al.l 
2006). The WMAP maps alone do not show any signif- 
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icant signature of the SZ effect (|Spergel et al.l [2006). Al- 
though the WMAP SZ signals of individual cluster have 
a low significance, it is possible to combine SZ signa- 
tures of many clusters to obtain constraints on the mean 
ICM properties. First efforts in this direction have been 
made through direct cross-correlation of galaxy /cluster 
surveys in optical/IR bands with the CMB tempera- 
ture maps obtained by WMAP dBennett et all 12 003a; 



Fosalba et al.ll2003l;lFosalba fc Gaztanagall2004l ; lMvers et all 



20041 ; I Afshordi et ai1l2004 ). These measurements have led 



to 2-5cr detection of anti-correlation with different galaxy 
and cluster surveys, which is consist ent with the expected 
SZ signal at WMAP frequencies dAfshordi et al.l |2004 



iHernandez-Monteagudo fc Rubino-Martfnll2004l ) 

More optimal detection of SZ effect is possible via cross- 
correlation of the CMB maps with the temperature of the 
X-ray emitting ICM, since both the SZ and X-ray probe the 
same hot ICM. This was a prim ary motivation for our study 
of the WMAP 1st year data. In lAfshordi. Lin fc Sanderson! 
i|2005l . ALS05 hereafter), we reported an 8a detection of the 
SZ signal, through an optimized correla tion of WMAP 1st 
year data release l|Bennett et al.ll2003bl ) with a catalog of 
X-ray clusters. Using a theoretically motivated ICM profile 
for the analysis, ALS05 have shown that the gas fraction 
in clusters is ~ 30 — 40% low compared to the universal 
baryon fraction of the Universe. If the clusters contain a 
repres entative mix of dark matter and bary o ns of the Uni- 
verse (lEvrardlll990l; IMetzler fc Evrardlll994l; INavarro et all 



Mohr et al 



1999; Bialck ct al 



2001 



19951; iLubin et al.lll99rj; lEke et al.lll99Sl; iFrenk et all 



Kravtsov et al.l 



1999 



2005 



Ettori et al.ll2006l : iMcCarthv et al.l200rj ). this number would 



be far too low to account for the expected baryon budgets 
in clusters, since the observed stell ar and cold gas frac- 
tion is no larger t han ~ 10 — 15% l|Gonzalez et al.l 120051 ; 
iLin fc Mohrl |2004). These results are in agreement with 
other SZ iLaRoque et~ai1 120061 ) and X-ray (|Ettoril 120031 ; 
IVikhlinin et al. 20061 ) studies. These results therefore indi- 
cate that the baryons are missing from the expected baryon 
budget of clusters. 

Much larger d iscrepancies with X -ray observations have 
been reported in iLieu et ah! l|2006l ). who claim that the 
WMAP SZ signal for ROSAT X-ray clusters is a factor 
of ~ 4 smaller than the expectation from X-ray bright- 
ness p rofiles. However, in order to compare with WMAP 
maps, ILieu et all ([2006) extrapolate the /?-model fits far 
beyond the region that is fit by X-ray data, suggest- 
ing that the reported discrepancy may simply be an ar- 
tifact of this e xtrapolation, which i s known to fail at 
large radii (e.g., IVikhlinin et al]|l999l ; ISchmidt et al. I l2004l ; 
IVikhlinin et al.ll2006l ). " 

Given the importance of the problem, we repeat the 
measurements of ALS05 using the 3 year WMAP data in 
this paper. The improvements from the previ ous analyses are 
two-f olds. First, the WMAP 3yr data release l|Hinshaw et al.l 
2006) should provide a higher significance detection of the 
SZ signal, while still capturing large angle information that 
is missing from higher resolution interferometric experi- 
ments. Second, we devise a model-independent method to 
measure the mean ICM pressure profile from WMAP tem- 
perature maps for a catalog of X-ray clusters. The latter will 
help minimize the dependence on theoretical assumptions 
and associated systematic uncertainties, which often domi- 
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Figure 1. The distribution of cluster redshifts and virial radii 
(estimated from X-ray temperature; see Sec. 12.21 1. The three upper 
lines show the resolution of WMAP bands (associated with the 
radius of t he disk with the sa me effective area as the detector 
beams; see Ijarosik et"aH l2006h . while the lower line shows the 
physical radius of the 1 degree circle at the cluster redshift. 



nate systematic uncertainties in SZ a nlayses carried out us - 
ing a cluster model (see e.g., ALS05: lLaRoaue et al.l l2006). 
We then use this measurement to put constraints on the 
thermal energy and baryonic budget of the ICM. 

The paper is organized as follows. In the next section 
we describe our methodology, including compilation of our 
X-ray cluster catalog, our method to extract the SZ signal 
from the CMB maps, and the hydrodynamic simulations 
that we compare to the data. Our main observational results 
are described in Sec. [3] In Sec. [4] we discuss the possible 
shortcomings, as well as the implications of our analysis, 
the most puzzling of which is a missing baryonic component 
of the ICM. Finally, Sec. [5] concludes the paper. 



2 METHODOLOGY 
2.1 X-ray cluster catalog 

Our primary objective is to study the SZ signal in a large 
sample of galaxy clusters by combining signals from known 
X-ray clusters. Therefore, we first assemble X-ray cluster 
catalog by searching the X-ray Galaxy Clusters Database 
(BAXJj. The database provides a comprehensive compi- 
lation of cluster properties from the literature, including 
the position on the sky, redshift, and X-ray measurements 
such as luminosity and/or temperature. The first pass in 
our cluster selection requires that clusters have measured 
X-ray temperature. We further remove clusters with mul- 
tiple entries in the database, and those that lie too close 



1 http://bax.ast.obs-mip.fr/ 
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to each other in projection. Although the resulting sample 
is not based on any statistically complete selection, we at- 
tempt to make the temperature measurements as consistent 
as possible, by adopting the ones that do not include the 
central regions of the clusters, or those obtained through 
a two-temperature component analysis, thus minimizing ef- 
fect of any "cool core" (CC) in the cluster mass determi- 
nation. Preference is given to a few studies that provide 
such temperature measu r ement for large sampl e s of clusters 
llFinoguenov et alJlioOll ; iFukazawa et alj|2004 llkebe et al.1 
2002). While in the former two studies, the central regions 
have been excised, in the latter, the temperature is given by 
the hotter component of their two-temperature model. This 
leaves us with a sample of 260 clusters with measured X-ray 
temperatures (Fig. [Q. 

For a fraction of our final cluster sample, estimates of 
the ICM gas cooling time i n the cluster centr a l region are 
available in the literatu re l|White et al.l 1 19971 : iPeres et al.l 
Il998l ; iBauer et al.ll2005h . thus allowing an investigation of 
the effect of the cool core on the cluster baryon budget. 
Cool cores are found in ~ 50% of nearby clusters, and al- 
though they show similar scaling properties across a wide 
mass range, they behave very differently from th e inner re- 
gions of non-CC clusters (|Sanderson et all l2006t ) . We refer 
to CC clusters as those with their central cooling times more 
than 3<r below half of the Hubble time (7 Gyr). In particular, 
for a subset of our sample with Tx > 3 keV (and after point 
source contamination cut, explained in Sec. l2,2|l . which com- 
prises 193 of our clusters, 80 clusters have measured central 
cooling times, and 39 of them are classified as cool core by 
this criterion. 



2.2 Extracting the SZ profile 

Our primary task is a measurement of the mean ICM pres- 
sure profiles by combining SZ signals from known X-ray clus- 
ters. In ALS05, we assumed a functional form for the ICM 
profiles to derive the pressure profile from data. In this work, 
we go a step further by developing a model-independent 
method to perform these measurements. This will help min- 
imize systematic uncertainties due to theoretical modeling 
of clusters, which are often the dominant sources of uncer- 
tainties in a model-dependent SZ analysis. 

The idea is based on the fact that clusters are ex- 
pected to b e self-similar and ex hibit universal dark mat- 
ter (NFW. iNavarro et al.l Il997h as well as ICM profiles 
l|Komatsu fc Seliakll200ll : lOstriker et alj|2005h . Recent cos- 
mological cluster simulations support that the ICM pres- 
sure profiles in cluster outskirt is remarkably self-similar, 
even in the presence of gas cooling and star for mation 
(|da Silva et al.l 120041 : iMotl et al.1 [20051 : iNagail 120061 ). Since 
the self-similar cluster forms a one-parameter family of the 
ICM pressure profile, both the radial extent and the ampli- 
tude of these profiles can be scaled by a single parameter, 
which can be chosen to be the total mass, X-ray or virial 
temperature of the cluster. 

The key ingredient of our method is to use the observed 
X-ray temperature of clusters to set the expected extent of 
the ICM profile, and use the WMAP maps to find the mean 
normalization. We use a simple analytic ICM model to set 
this scale: 



*»° = ( L16 MPC > ( tOOkm/s/Mpc ) 1 (sSv) V2 ' (1) 

where 7-200 is the radius of the sphere within which the mean 
density of the cluster is 200 times the critical density of the 
Universe, p cr it = 3H 2 (z) / (8irG) , H(z) is the Hubble con- 
stant at the redshift of the cluster, and Tx is the observed X- 
ray temperature of the cluster. Following ALS05, the model 
is based on a spherically symmetric NFW gravitational po- 
tential with concentration C200 = 5, and a polytropic ICM in 
hydrostatic equilibrium. Note that, for a cluster of a given 
Tx, Eq. [T] provides r2oo that agrees with the best fit scaling 
relation of the recent Chandra observation of X-ray clusters 
l|Vikhlinin et al.l 120061 ) with the accuracy better than 5% 
leveE 

Once the profiles are scaled at r2oo, we search for the 
SZ signal out to 4r2oo radius, since the shock heated gas is 
only expected to be present within a few times the virial ra- 
dius of the clusters. Given that the CMB signal is correlated 
on the degree scale (Doppler peak) in the sky, it is impor- 
tant to distinguish between primary CMB and SZ signal, 
and also account for WMAP beam smearing. Therefore, we 
include pixels out to 8^200 + 2 deg, where #200 is the angular 
size of r2oo in the skj0. We then aim to constrain spheri- 
cally averaged pressure profiles within radial bins that are 
logarithmically spaced in radius, and are centered at 0.25, 
0.5, 1, 2, and 4r2oo- We assume that pressure is smoothly 
interpolated as P ga s(r) = A + Br~ 3 in between the centers 
of the bins, where A and B are constants, and r~ 3 behav- 
ior is motivated by the dark matter/gas/pressure profile in 
the outskirts of simulated haloes (see Fig. [2}. For complete- 
ness, we further assume P ga , B (r) oc r~ 2 within 0.25r2oo of the 
cluster center. 

The dominant sources of noise for our SZ measure- 
ment on angles that can be resolved by WMAP experiment 
are WMAP detector noise and CMB primary anisotropies, 
which are both expected to be Gaussian to a good approxi- 
mation. Therefore, minimizing 

X = 51 [Tia - Sia] C^.jb [Tjb - Sjb] , (2) 

is the optimum method of constraining the cluster SZ pro- 
file, Sia (which is a superposition of our radial pressure 
bins), from the observed CMB temperature maps T a . Here 
i and j sum over WMAP pixels, while a and b sum over 
WMAP frequency bands, Q(41 GHz), V(61 GHz), W(94 
GHz). The noise covariance matrix d a ,jb is the sum of pri- 
mary CMB correlation function (using CAME0; Lewis et al.l 
2000), which dominates on large pixel separations, and the 
WMAP detector noise which is assumed to be uncorrelated 



2 However. iNaeai et al. I (1200711 find that the estimated r2oo from 
hydrostatic equilibrium is biased low (by about 10%) on average, 
if subsonic turbulent gas motions are present and t heir pressure 
contr ibution is not explicitly accounted for (also see iRasia et al.l 
2006). In Sec. 14.21 we discuss the implications of using a higher 
normalization for the r200 — Tx relation. 

3 For example, for a small cluster, in order to find the SZ profile, 
the mean (large angle) primary CMB background needs to be 
subtracted. Including pixels beyond the cluster helps with setting 
this background. 

4 |http://camb.info/| 
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for different pixels/frequencies. Note that both Ci a .jb and 
Sia need to be convolved with detector beam+pixel window 
functions. For WMAP, the effective beam radii range from 
0.3 deg for Q-band to 0.1 deg for W band, while we use 
WMAP foreground cleane d maps in N S ide = 512 HEALPix 
format (|G6rski et al.|[2005l ; (0.1 deg) 2 pixels). We mask out 
all the pixels within Kp2 Galactic mask (~ 13% of the sky), 
but do not mask out the point sources, as they are corre- 
lated with the clusters in our sample. In order to deal with 
a tractable covariance matrix, we progressively degrade the 
map pixel resolution towards the cluster outskirts, so that 
we have < 800 pixels per cluster. Moreover, in order to ac- 
count for the uncertainty in measured Tx, as well as a 15% 
intrinsic scatter in the r2oo — Tx relation (seen in our simu- 
lated clusters; see Sec. l2.3p . we divide the x 2 for each cluster 
by a constant factor. 

Point sources are additional sources of contamination 
of the SZ signal. The contamination is expected to be the 
largest in the lowest frequencies. So, to minimize the contam- 
ination, we first exclude any cluster that contains a radial 
bin with larger than +3a Q-band signaQ (remember that 
SZ signal is negative for WMAP frequencies). This excludes 
18 of our clusters, only 5 of which have Tx > 3 keV. Ad- 
ditionally, we assume a radio point source with a fiat spec- 
trum, which is shown to be consist ent with the spectrum of 
point sources observed by WMAP (jBennett et al.ll2003al ) , at 
the center of each cluster. Note, however, that our results 
remain virtually unchanged even if we use a steeper spec- 
trum obtained in lCoble et al.l (|2006h for fainter cluster radio 
sources. To find the ICM pressure profile, we then marginal- 
ize over the mean amplitude of the central source, as well 
as a mean point source contamination of the innermost ra- 
dial bin, which are expected to be the dominant sources of 
contamination of the SZ signal. 
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Figure 2. Mean pressure profile of 193 of our most massive clus- 
ters with Tx > 3 keV (points+errorbars). The gray area shows 
the 68% region allowed by this measurement, assuming P ga s > 
prior, which reflect the errors, as well as their correlations. The 
blue/solid curves are predicted pressure profiles from nine simu- 
lated clusters with Tx > 3 keV, while black/dotted curves show 
dark matter density from the same simulations, divided by the 
critical density of the Universe. 



2.3 Hydrodynamic Simulations 

In this study, we use high-resolution cosmological simula- 
tions of nine hot {Tx > 3 keV), massive galaxy clusters 
forming in the flat ACDM model: fi m = 1 - tt A = 0.3, 
fib = 0.04286, h = 0.7 and a 8 = 0.9, where the Hubble 
constant is defined as lOO/i km s _1 Mpc -1 , and erg is the 
power spectrum normalization on 8/i -1 Mpc scale. The sim- 
ulations follow dissipationless dark matter and stars, as well 
as dissipative gas dynamics, and include a number of phys- 
ical processes critical to galaxy formation, including radia- 
tive cooling, star formation, stellar feedback and metal en- 
richment. The simulations were performed with the Adap- 
tive Refinement Tree (ART) jV-bod y+gas dynamics code 
l|Kravtsovll999l ; [Kravtsov et alj2002h . an Eulerian code that 
uses adaptive mesh refinement to greatly increase the reso- 
lution in the high density region and resolve formation and 
evolution of cluster galaxies and their impact on cluster gas. 
Massive clusters with Tx > 3.7keV, for example, are simu- 
lated in a box size of 120 Mpc with the spatial resolution 
of 3.5/i _1 kpc and the mass resolution of 9 x 10 8 /i _1 Mq. 
Less massive clusters are simulated in a smaller box of 
80 Mpc in size with higher spatial (~ 2.5 /i~ 1 kpc) and mass 

5 Note that the systematic bias in the mean due to a 3<r trunca- 
tion of a Gaussian probability function is less than 0.005cr. 



(3 x 10 8 /i -1 Mq) resolutions. Throughout this work, we use 
two different average temperatures of simulated clusters. 
The first is a X-ray spectral temperature, Tx, a value de- 
rived from a single-temperature fit to the integrated clus- 
ter spectrum within rsoo extracted from the mock Chan- 
dra images excludin g the core (r < 0.1 5r5oo) and detectable 
small-scale clumps (Na gai et alj|2007h . We also compute a 
gas-mass-weighted temperature, denoted as (T) throughout 
this paper, obtained by weighting of the 3D temperature 
with the gas density measured directly in simulations. Note 
that (T)200 is the gas-mass- weighted average temperature 
measured within r^oo- These simulations have been used 
to study the effects of galaxy formation on the Sunyaev- 
Zel'dovic h effect ([N agai 2006) and the baryon fraction in 
clusters l|Kravtsov et al.ll2005l ). These studies are based on 
a sample of 11 simulated clusters with only 4 clusters mas- 
sive enough to be included in the current study. Since then, 
we have simulated five more massive clusters, making the 
total sample of 16 clusters, from which we select nine mas- 
sive clusters with Tx > 3 keV for the current study. The 
detailed descriptio n and our cluster sample are described in 
l|Nagai et alj|2007t ). 

Note that, in order to compare simulations with obser- 
vations, the simulated gas density or pressure is normalized 
by the ratio of cosmic con cordance baryonic m ass fraction 
(= D. b /Q m = 0. 168 ± 0.007; ISpergel et~aLll2006l ) to the value 
used in the simulations. 
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Figure 3. The ratio of total thermal energy to the mass times 
Tx , for the same clusters as in Fig. [2] The colors are the same as 
in Fig. [2] and we integrate over the pressure profiles to find total 
energy, M gas (T)(< r), enclosed within radius r, while we esti- 
mate the total mass from an NFW profile. The dot ted lines show 
the t otal cosmic baryonic budget of 0.168 ± 0.007 jSpergel et al.l 
2006). The dashed region show s the Chandra observa tional con- 
straints from 8 X-ray clusters l|Vikhlinin e t al. 2006j), while the 
open point+errorbar is the most recent SZ con straint from 
OVRO/BIMA interferometers l lLaRoque et alj|2006h . 



3 OBSERVATIONAL RESULTS 
3.1 Universal Pressure Profile 

The main result of our analysis is shown in Fig. [2] It shows 
the mean ICM pressure profile (P = p gas T), normalized 
to the critical density of the Universe at the cluster's red- 
shift times its observed X-ray temperature (p CI itTx)- In 
this plot, we only include clusters with Tx > 3 keV (or 
A/200 > 2.4 x 10 14 M Q ). The points+errorbars show our best 
fit measurement for a universal pressure profile. The gray 
area shows the 68% region allowed by this measurement, 
assuming P gas > prior, which reflect the errors, as well 
as their correlations. The best fit measurement of Fig. [2] is 
preferred to null at Ay 2 = 47. The simulated profiles can 
improve % 2 by 31 — 43, which implies a 6 — 7a confirmation 
of the simulated models, with respect to null. 

Our measurements of the mean ICM pressure profiles 
closely match the numerical prediction for the pressure pro- 
files of nine simulated clusters (indicated by blue/solid lines) 
in the same temperature range. It is also interesting to notice 
the similarity of the observed ICM pressure and the simu- 
lated dark matter density profiles, indicated by black/dotted 
lines. 



3.2 Total Thermal Energy of the Cluster 

In Fig. ©, we examine the integrated pressure, or equiva- 
lently themal energy of the ICM as a function of cluster cen- 
tric radius in units of r2oo- Here we normalize the thermal 
energy (M gas ,2oo (P)20o) to the total mass of the cluster (ex- 



pected from an NFW profile with C200 = 5) times the cluster 
X-ray temperature (M tot ,20oTx). The solid point+errorbar 
shows the mean and standard deviation of this quantity at 
T200, and thus corresponds to our measurement of the mean 
ICM thermal energy: 



gas, 200 



tot, 200 



<r> 2 



= 0.095 ±0.011. 



(3) 



The results is in very good agreement with our simulated 
clusters 



M, 



gas, 200 



Mt 



tot, 200 



<r> 2 



Tx 



= 0.093 ±0.003, 



(4) 



which are shown by thick curves in Fig. [3] Finally, the dashed 
region sho ws the most detailed available X-ray study of 
8 clusters (|Vikhlinin et al.l 120061 ) (T x > 3 keV), observed 
by Chandra X-ray observatory, while the open dot with 
the errorbar shows the constraints from SZ observations 
of 38 massive cluster s with OVRO/BIMA interferometers 
l|LaRoaue et al.ll2006T ). Both of these studies are also fully 
consistent with our SZ measurements from the WMAP 3 
year data. 

Knowing (T)20o/2jf, we can convert the thermal energy 
fraction into the gas fraction of clusters. For our nine simu- 
lated clusters with Tx > 3 keV, we estimate that the ratio 
of mean to X-ray estimated temperature is (T)2oo/Tx = 
0.87 ± 0.02. Using this value, the cluster gas fraction is 



S.200 



M t< 



0.109 ±0.013 



(•») 



for the WMAP data. The value is significantly smaller 
than the cosm i c bar yon fraction of f2f,/f2 m = 0.168 ± 0.007 
|Spergel et alj [2006). If the clusters contain a represen- 
tative mix of dark matter and baryons of the Universe 
(Mbaryons/Mtot ~ fi&/f2 m ), our results indicate that a frac- 
tion 



Af, 
Mi 



111 1 1 11 \ 



M s 



baryons 



Mi 



= 0.35 ±0.08 



baryons 



(6) 



of the cluster baryons are missing from the hot ICM. As we 
discuss more in § 14.31 this number is far too large to be ac- 
counted for by the obser ved stellar and cold gas fraction of 
« 10 % in our T x range (|Gonzalez et alj|2005t iLin fc Morirl 
2004). These results are in agreement with our earlier anal- 
ysis by ALS05. 



3.3 Systematic Trends with Tx 

Exploiting a large coverage of cluster X-ray temperatures 
and masses, it would be interesting to study if there is any 
systematic trend in the gas fraction as a function of Tx- 
Fig. [4] shows the thermal energy to mass ratio within r2oo 
as a function Tx , where we have binned clusters into 2 keV 
temperature bins. The number of clusters in each bin, from 
left to right is 44, 66, 84, 41, 16, and 9. Here again, we see 
that our SZ measurements are consistent with predictions 
from hydrodynamic simulations (open stars), as well as our 



6 The quoted error is the sample variance error of the mean en- 
ergy fraction for the simulated clusters with Tx > 3 keV. 
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Figure 4. Thermal energy to mass times Tx ratio as a func- 
tion of Tx • Solid points+errorbars show the measured total ther- 
mal energy to mass ratio within r2oo for all our clusters, binned 
within ATx = 2 keV bins. Open stars show the same ratio for 
our hydrodynamic simulations, while the gray area shows the 68% 
uncertainty range of the best power law fit to our data (Eq. 0. 
The black h orizontal bar shows the total cosmic baryonic budget 
l|Spergel et al . 2006). multiplied by the mean tempcratr 
tion from the simulations. 



turc coitcc- 



earlier analysis (ALS05). The best fit power law to our mea- 
surement points (weighted by their errors) is 



gas, 200 



Ml 



tot, 200 



T x 



(0.098±0 



015) ( 



6.4 keV 



which is represented by the gray area in Fig. [4] and does 
not show any significant monotonic trend throughout the 
Tx range, a trend similar to the simulated clusters. 

The measurement points in Fig. [4] have almost a factor 
of two more scatter than is expected from the measurement 
errors. This could be caused by a possible systematic un- 
derestimate of our measurement errors, or otherwise by an 
intrinsic scatter in the ICM properties of different clusters. 
An important clue to the reason behind this scatter, as well 
as the nature of missing baryons may come from looking at 
a sub-sample of our clusters with resolved cool cores (39 CC- 
clusters with Tx > 3 keV; see Sec. 12. ip . The energy fraction 
of the CC-clusters is 



M s 



M u 



{Th 



T x 



= 0.142 ±0.019, 



(8) 



which is significantly larger than the main sample, and 
is consistent with no missing baryons after the correction 
for {T)2oo/Tx ■ Moreover, the pressure profile seems signifi- 
cantly flatter than the simulated profiles within 0.5r2oo- Al- 
though the fact that the ICM gas fraction is systematically 
higher for cool core clusters in Eq.[8]may signal an unknown 
systematic problem with our method of subtracting point 
sources, it may also be caused by a real anti-correlation be- 
tween the missing baryons and the central cooling activity. 
Further observational studies are clearly needed to resolve 
these issues. 



Although our sample includes X-ray clusters with mea- 
sured temperatures out to redshift 1, less than 20% of the 
statistical significance of our SZ detection in our sample 
comes from clusters with z > 0.2. This is because the virial 
radius of clusters (~ 1.5 Mpc) at z > 0.2 falls below the 
WMAP detector beam size of ~ 0.2 deg. Therefore, it would 
be difficult to study any redshift evolution in our sample, 
and thus it will not be explored here. 



4 DISCUSSIONS 

4.1 Comparison with our previous analysis 

One may wonder why, despite the use of the 3 year WMAP 
maps as well as more X-ray clusters, the significance of our 
SZ detection has not significantly increased compared to our 
previous analysis with the 1st year WMAP maps in ALS05. 
This is mainly due to the fact that, in our new analysis we 
parameterize the mean cluster signal by 8 free parameters (6 
for pressure profile, and 2 for point sources), while our pre- 
vious analysis only used 2(1 for the total gas fraction, and 
1 for the over all point source contamination). Therefore, 
the error on the individual parameters, such as the total 
gas fraction, is increased due to the extra marginalizations 
involved in the analysis. 

Despite the use of a completely different method, our 
main conclusion: a fraction of 30% — 40% missing ICM 
baryons, has not changed from ALS05. However, we did not 
find any systematic Tx dependence in the energy fraction in 
Eq. [7] This shows a ~ 2<r discrepancy with our earlier anal- 
ysis of WMAP SZ signal, where we found an almost linear 
dependence on Tx ■ However, our earlier analysis assumed a 
universal pressure profile for all clusters. Therefore, one way 
to reconcile the two analyses is for the profiles to be sys- 
tematically shallower for clusters with Tx ^ 4 keV, which 
is indeed seen in our sample, but it is difficult to quantify 
due to the low significance of SZ signal in cooler clusters. 



4.2 Systematic and Theoretical Uncertainties 

One of the main improvements from our previous analy- 
sis in ALS05 is the use of a model-independent method to 
derive the mean ICM profiles from the WMAP 3yr data. 
This will help minimize systematic errors arising from un- 
certainties in theoretical modeling, which often dominate 
systematic uncertainties of the analyses carried out using a 
model-dependent method. 

The only theoretical assumption in our derivation of 
the mean ICM pressure profile has been the use of Eq. [T] 
to infer the cluster masses or radii based on their mea- 
sured X-ray temperatures. Therefore, any systematic errors 
in this scaling relation translate into an error in the cluster 
total mass estimate, which could bias our estimate of the 
cluster gas mass fractions. Indeed, the normalization of the 
7*200 — Tx used in our analysis is con sistent with hydrostati c 
estimates from X-ray observations l|Vikhlinin et al.l 120061 ) , 
which a re about 10% smalle r than the simulated normal- 
ization |Kravtsov et alj [20061. However, even if there is a 
systematic error in the normalization of r2oo — Tx relation, 
Eq.[l]can be viewed as the definition of r2oo in terms of the 
observed cluster X-ray temperature, which simply defines a 
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unified radial unit for the mean pressure profile. In other 
words, the bias, if present, results in the re-definition (or 
re-scaling) of r2oo- In this way, our measurements of mean 
pressure profile is indeed model-independent. However, any 
inference about the missing ICM baryonic fraction would be 
affected by a systematic error in our assumed r2oo — Tx rela- 
tion. For example, using the higher normalization (inferred 
from simulations) can boost the missing baryon fraction to 
~ 50%. 



4.3 Missing Baryons or New Astrophysics? 

The most puzzling result of our analysis is that a significant 
fraction (35% ± 8%) of the baryons are missing from the 
ICM. On average, the fraction of the observed stellar mass 
associated with galaxies in clus ters within our Tx range is 
less than 10% (|Lin et all [20031 ). leaving room for an un- 
accounted baryonic component of 25 ±8%. Although the ev- 
idence for this component is marginal (3.1er) from our analy- 
sis, inde pendent verification from other SZ a nd X-ray obser- 
vations (|LaRoaue et all 120061 ; lEttoril 120031 ; IVikhlinin et al.l 
120061 . see Fig. |3]) strongly support the case for a missing 
component of cluster baryons. 

Interestingly, modern cosmological cluster simulations 
that include gas cooling and star formation also give the 
ICM mass fraction v ery similar to the observed values 
l|Kravtsov et all [20071) . The good agreement in the ICM 
mass fraction between simulations and observations suggests 
that the simulations reproduce the hot ICM of the real clus- 
ters reasonably well ( Nagai et al.ll2007h . However, the agree- 
ment is achieved while the stellar fractions within r^oo of the 
simulated clusters are a factor of two to three greater than 
the observed stellar fraction. 

This raises a possibility that a significant fraction of 
baryonic mass may be hidden in a form that is difficult 
to detect with ordinary observational means. One such ex- 
ample is diffuse intracluster light (ICL), which would be 
challenging to detect observationally, but could potentially 
hide significant amount of stellar mass in the vast intra- 
cluster space. The measurements of ICL are now becom- 
ing possible with very deep photometry. Although the mea- 
surements are quite challenging, a number of ICL observa- 
tions indicate that no mor e than 10-50% of stars could be 

Lin fc Mohrll2004 IZibetti et all 



Monaco et al.ll2006l ). Similar re- 



found in_theJbrm_of ICL 
l2005l : lGonzalez et al.l[2005l : 

suits are also obtained in cosmological cluster s imulations 
|Murante et al.ll2004l ; ISommer-Larsen et al.ll2005l ). Alterna- 
tively, one may also envisage hiding baryons in a form of cold 
and compact dark baryonic clouds formed from local cooling 
instabilities within ICM (or proto-ICM), similar to the high 
velocity clouds in the local group. At the moment, there 
is essentially no constraint on such clouds as long as their 
temperature is below 10 K, and/or their covering fraction of 
the c luster radio map is less than unity (Dwarak anath et al.l 
1 19941 ). 

It has been suggested that a soft X-ray excess in obser- 
vations of X-ray clusters can be interpreted as evidence for a 
warm (~ 0.1 keV) compact phase of the ICM, which would 
not contribute to the SZ signal, while co ntaining a signifi- 
cant fraction of the ICM baryons (e.g., iBonamente et aD 
120051 . and references therein). However, the extragalac- 
tic nature of the soft X-ray excess is still controversial 



(e.g., iBregman fc Llovd-Daviesll2006l ). Moreover, it is not 
clear how the strong cooling instability of any such compo- 
nent could be avoided in the cluster environment. 

Yet another pos sibility that has been recently pro- 
posed by lLoebl (|2007f ), is thermal diffusion or evaporation 
of baryons out of the virial radius of the cluster. The ther- 
mal evaporation, if efficient, can remove the gas out of the 
virial radius, bringing the cluster gas fracti on down with- 
out creating stars (but see lMedvedevl [2007). Similarly, in- 
clusion of diffusive processes in hydro dynamical simulations 
may on one hand suppress cooling (jZakamska fc Naravanl 
l2003h . thus bringing down the cool gas fraction into better 
agreement with the observed stellar content, while on the 
other hand r emove the t hermal energy of the ICM out of the 
virial radius (|Loebll2002t ). explaining the low gas fractions in- 
ferred from our SZ observations. However, recent numerical 
simulations that included (some of) these processes suggest 
that they cannot a ffect the overcooling problem significantly 
l|Dolag et al.ll2004l ). 



5 SUMMARY AND FUTURE PROSPECTS 

We report a 9cr detection of the SZ signal in the CMB maps 
of WMAP 3yr data, through study of a sample of 193 mas- 
sive clusters with observed X-ray temperatures larger than 
3 keV. The improvements from the previous analyses in 
ALS05 are in both data and ana lysis methods: (1) th e use 
of the WMAP 3yr data release (fflnshaw et al-lboOfJ ) and 
(2) the use of a model-independent method to measure the 
mean ICM pressure profile to minimize systematic uncer- 
tainties due to theoretical modeling. 

The resulting pressure profile is in good agreement with 
measurements of recent X-ray and SZ observations as well 
as hydrodynamical cluster simulations. Our result indicates 
that 35 ± 8% of the baryons are missing from the ICM, 
which is significantly larger than the observed stellar or cold 
gas fraction of about 10% so far accounted for in clusters. 
The evidence for this component is marginal (3.1ct) from 
our analysis, but ind ependent verification fr om other SZ 
llLaRoaue et al.ll2006t ) and X-ray observations (|Ettorill200"3l : 
IVikhlinin et al.ll2006l ) strongly support the case for a missing 
component of cluster baryons. This signals the presence of a 
missing baryonic component or yet-unknown astrophysical 
processes that could lower the cluster baryon fraction. 

The problem of missing thermal energy of the ICM in- 
troduces a new puzzle into the standard cluster formation 
scenario. Further investigations of the ICM physics are thus 
essential for our understanding of cluster formation as well 
as applications of cluster-based cosmological tests using X- 
ray and SZ cluster surveys. On the theoretical front, it is im- 
portant to understand details and relative importance of gas 
cooling and heating by supernovae (SN) and active galac- 
tic nuclei (AGN) and their effects on the ICM properties. 
Further numerical investigations of diffusive processes are 
also critical for understanding of cluster plasma phenomena. 
Observationally, it is important to further refine measure- 
ments of X-ray, SZ, lensing and optical/IR observations and 
consolidate the measurements of the cluster baryon budgets 
in clusters. Large microwave interferometer arrays such as 
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will ultimately be able to map the ICM pressure 
profile in the outskirts of ind ividual clusters with exquisite 
accuracy (jKocsis et al.ll2005h . while the next generation of 
high-resolution, multi-frequency CMB experiments, such as 
APEX, ACT, SPT, and SZA0, will be able to probe the ther- 
mal energy content to high redshifts, opening new windows 
for the future studies of structure formation and cosmology. 
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